The Clustering of AGN and Galaxies at Intermediate Redshift 
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ABSTRACT 

Galaxies in the environments of 69 0.2 < z < 0.7 UBR selected AGN have been imaged 
to Bj ~ 23.5. By applying photometric redshifts and color selection criteria to the galaxy 
catalogue, the AGN-galaxy cross-correlation function has been measured as a function of galaxy 
type. The spatial cross-correlation of AGN with red (early-type) galaxies is comparable to the 
autocorrelation function of elliptical galaxies at low redshift. In contrast, the cross-correlation 
of AGN with blue (late-type) galaxies is weak and has been detected with low significance. As 
blue galaxies dominate Bj ~ 23.5 galaxy catalogues, the cross-correlation of UBR selected AGN 
with all galaxies is weak at intermediate redshifts. 
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The association of Active Galactic Nuclei hosts 
with other galaxies is well established (Bahcall et 
al. 1969; Yee and Green 1987; HaU and Green 
1998; Laurikaninen and Salo 1995; Smith et al. 
1995, 2000) and provides constraints on the mod- 
els for the formation and fueling of AGN. In ad- 
dition, as QSOs can be used to trace large-scale- 
structure at z 1 (Boyle et al. 1999), estimates 
of the QSO environments are required to estimate 
their bias with respect to galaxies which are used 
to trace large-scale-structure at z < 1. 

Several previous studies of AGN environments 
are summarized in Tabic 1. Radio- loud AGN 
are typically found in environments comparable 
to galaxy clusters (Yee and Green 1987; Hall and 



Green 1998; Wold et al 2000) while z < 0.3 radio- 
quiet AGN appear to be associated with poorer 
environments similar to field galaxies (Laurikani- 
nen and Salo 1995; Smith et al. 1995; De Robertis 
et al. 1998). While the z > 0.3 radio-quiet AGN- 
galaxy cross-correlation function is comparable to 
the clustering of field galaxies (Smith et al. 2000), 
it has only been detected with 2a significance 
(EUingson et al. 1991; Teplitz et al. 1999; Smith et 
al. 2000) and no detailed information is available 
on the strength or evolution of galaxy clustering 
around radio-quiet AGN at these redshifts. 

Previous studies of AGN environments have 
typically used galaxy catalogues derived from 
single band imaging. Galaxy evolution and 
fc— corrections result in a changing morphologi- 
cal mix of galaxies as a function of redshift and 
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depth. As galaxy clustering is eorrelated with 
morphology and color, the clustering properties 
of galaxies are also a function of redshift and the 
observer's bandpass. It is therefore plausible that 
estimates of AGN environments at z > 0.3 have 
been biased by the changing properties of the 
galaxy catalogues selected from single band imag- 
ing. By using color selection, it is possible to select 
the same galaxy type as a function of redshift. By 
measuring the cross-correlation function of AGN 
with early and late-type galaxies, it is possible 
to determine if radio-quiet AGN are in unusual 
environments. 

2. Data 

2.1. The Galaxy Sample 

The galaxy sample has been previously used by 
Brown et al. (2000) to measure the clustering of 
galaxies as a function of color and a more detailed 
description of the catalogue is provided by Brown 
et al. (2000) and Brown (2000). The image data 
consists of 5° x 5° images of the South Galactic 
Pole (SGP) and UK Schmidt field 855 (F855). The 
images were produced by stacking SuperCOSMOS 
scans of UK Schmidt photographic plates in U, 
Bj, R and I bands. Object detection, instrumen- 
tal photometry and faint object star-galaxy classi- 
fications were determined with SExtractor Bertin 
and Arnouts (1996). Photometric calibration of 
the data was determined with CCD images and 
published photometry. To prevent dust extinction 
introducing spurious large-scale-structure, magni- 
tude estimates are corrected with the extinction 
estimates of Schlegel et al. (1998). The final 
galaxy catalogues are complete to Bj ^ 23.5 and 
contain ~ 2 x 10^ galaxies per field. 

2.2. The AGN Sample 

The AGN sample consists of 69 0.2 < ^ < 0.7 
UBR selected broad emission line AGN from La 
Franca et al. (1999). The survey area is the SGP 
field and while the catalogue does not have homo- 
gcneoiis sky coverage, it is not strongly concen- 
trated in any one part of the field. As the F855 
field only contains 7 Q.2 < z < 0.7 AGN with pub- 
lished positions (Vcron-Cctty and Veron 2000), 
the F855 field has not been used to measure AGN- 
galaxy clustering. AGN positions were determined 
by selecting Bj < 23.5 objects in the stacked 



SuperCOSMOS scans within 5" of the published 
positions. The resulting catalogue of AGN con- 
tained 69 Bj < 21 objects with U - Bj < 0.1 
colors. The Bj absolute magnitudes and redshifts 
of the AGN are plotted in Figure 1. The cata- 
logue contains 36 QSOs (Mbj < -21.5 + 51og/i 
where h x 100kms"^Mpc"^ = Hq) and 33 Scyfert 
1 galaxies. Five of the AGN are within 20" of 
1.4 GHz sources detected by the NVSS (Condon 
et al. 1998) but all have fluxes less than 60 mJy 
and 1.4 GHz luminosities less than lO^^WHz"^. 

2.3. Photometric Redshifts 

Approximately 700 < z < 0.8 galaxy red- 
shifts in the SGP and F855 fields are available 
from the NED database. Spectroscopic redshifts 
are available for galaxies as faint as Bj ^ 24 
so it possible to use multicolor photometry and 
spectroscopic redshifts of galaxies in the SGP and 
F855 fields to calibrate the photometric redshifts. 
The relationship between the multicolor photom- 
etry and the galaxy redshift was determined by 
fitting quadratic functions to the data (Connolly 
et al. 1995). The relationship was determined for 
UBjRI, UBjR, BjRI and BjR photometry as 
only a small fraction of the catalogue is detected 
in all 4 bands. Figure 2 shows a comparison photo- 
metric redshifts and the spectroscopic redshifts for 
galaxies in the two fields. While BjR photomet- 
ric redshifts are poorer than photometric redshifts 
derived in 3 or more bands, they do place useful 
constraints on the redshifts. Error estimates for 
the photometric redshifts have been determined 
by the measuring the rms of the residuals as a 
function of photometric redshift and color. The 
accuracy of the photometric redshifts is a func- 
tion of galaxy color and, as shown in Figure 3, red 
galaxies have comparatively small errors. 

2.4. Color Selection 

A significant bias present in most studies of the 
AGN environment is that they have relied on small 
samples or single band imaging. Galaxy evolution 
and fc— corrections result in a changing morpho- 
logical mix of galaxies as a function of limiting 
magnitude with single band imaging. Galaxy cat- 
alogues in bands bluer than R are dominated by 
weakly clustered (blue) late-type galaxies at mag- 
nitudes fainter than Bj ~ 22 (Efstathiou et al. 
1991). 
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Multicolor imaging provides significant advan- 
tages for the study of the AGN host environment. 
Colors and photometric redshifts can be used to 
select particular galaxy types at faint magnitudes. 
As galaxy clustering is strongly correlated with 
color (Brown et al. 2000) and morphology (Davis 
and GcUcr 1976; Loveday et al. 1995), studying the 
cross-correlation of AGN with particular galaxy 
types should help determine if AGN hosts are in 
unusual galaxy environments. 

The color criteria applied to the galaxy cata- 
logue select galaxies redder or bluer than a non- 
evolving Sbc in Bj — R where the color as a func- 
tion of redshift is determined with fc— corrections 
from Coleman et al. (1980). The use of Bj - R 
rather than shallower U — Bj or R — I allows 
the selection oi Bj = 23.5 early- type galaxies at 
~ 0.5. As the Bj—R selection criteria are a func- 
tion of redshift, photometric redshifts are used to 
determine the correct value of Bj — R when se- 
lecting galaxies. For the remainder of this paper, 
the early and late subsaniplcs will refer to galax- 
ies redder and bluer than the Bj — R selection 
criteria. Figure 4, a plot of the colors and photo- 
metric redshifts of galaxies in the SGP with HST 
morphologies (Abraham et al. 1996; Small et al. 
1997), shows that the color selection criteria are 
capable of selecting different types of galaxies. 

3. The Angular and Spatial Correlation 
Functions 

The angular correlation function estimates the 
fractional excess of object pairs at a given angu- 
lar separation compared with what would be ex- 
pected for a random distribution of objects. The 
estimators of the angular correlation function used 
in this paper require random object catalogues 
which are produced by making copies of the galaxy 
catalogue and randomly repositioning the copied 
galaxies across the field. 

The angular autocorrelation function of galax- 
ies is estimated with 

DD-2DR + RR 
= RR 

(Landy and Szalay 1993) where DD, DR and RR 
are the number of galaxy-galaxy, galaxy-random 
and random-random pairs at angular separation 
6 ± 66. The Landy and Szalay (1993) estimator is 
only applicable to the autocorrelation function so 



the AGN-galaxy (and early-late) cross-correlation 
function is determined with 



(2) 



where AG and AR are the number of AGN- 
galaxy and AGN-random pairs at angular sepa- 
ration 8 ±69. The random errors of AR, DR and 
RR are reduced by taking the average of multi- 
ple estimates of each parameter where each esti- 
mate has been determined with a different random 
object catalogue. Both estimators of the angular 
correlation function satisfy the integral constraint. 
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(Groth and Peebles 1977), resulting in an midcr- 
estimate of the angular correlation function. To 
remove this bias from the correlation function, the 
term 

io{e)n = -^J J ^{S)S^iS^2 (4) 

is added to the estimate of the correlation func- 
tion. The term u{9)q does require an assumption 
of the form of the correlation function to correctly 
estimate the value of correlation function. 

Galaxy angular correlation functions are typi- 
cally approximated by power-laws where 



^(6*) = Ae^-'^ 



(5) 



where ^4 is a constant and 7 ^ 1.7. The corre- 
sponding spatial correlation function is given by 



^{r,z) = {r/ro)-^l + z)-^^+'^ 



(6) 



where r is the separation of the galaxies in physical 
coordinates and ro and e are constants. For e = 
and e = 7 — 3, the clustering properties are fixed 
in physical and comoving coordinates respectively. 
The relationship between the angular and spatial 
correlation functions is given by Limber's equation 
(Limber 1954). If ^(r, z) ^ when r > O.lz, then 
Limber's equation is given by 



(J 



/I'tdzJ^dz 



dz 



(7) 



(PhiUipps et al. 1978) where ^ and ^ are the 
redshift distributions of each set of objects (e.g. 



3 



the early and late subsamples) and r{0) is the dis- 
tance between two objects at z and z' separated 
by angle 6 on the sky. For the AGN-galaxy cor- 
relation function, the rcdshifts of each AGN are 
known so Limber's equation can be written as 



^{0)= Er HaMO),z)'-^dz 



dz 



dz 



(8) 



where 71 is the number of AGN, is the num- 
ber of galaxies per unit redshift and r{6) is the 
distance in physical coordinates between an AGN 
at redshift Zi and a galaxy at redshift z separated 
by angle 9 on the sky. 

4. The Clustering of Faint GalcLxies 

Estimates of the clustering of z ^ 0.5 galaxies 
are required if the environment of AGN hosts is 
to be compared to the "normal" galaxy environ- 
ment. To make the comparison valid, the same 
galaxy selection criteria are applied to the study 
of galaxy-galaxy clustering as AGN-galaxy clus- 
tering. The early and late subsample autocorrela- 
tion functions and the early-late cross-correlation 
functions have been determined in the SGP and 
F855 fields. Plots of the 18.0 < Bj < 22.5 an- 
gular correlation functions are shown in Figure 5. 
The amplitudes and values of the 7 as a function 
of magnitude are summarized in Tables 2 to 4. 
The value of 7 and the amplitude of the cluster- 
ing strongly depend on color (Brown et al. 2000). 
The observed clustering in the SGP is consistently 
stronger than the observed clustering in F855 at 
bright magnitudes. However, at fainter magni- 
tudes the difference between the clustering prop- 
erties in the two fields is significantly reduced and 
the late subsample autocorrelation functions are 
comparable. At faint magnitudes, the estimates 
of 7 for the early-late cross-correlation function in 
the SGP and F855 fields differ by ~ 2a. It is pos- 
sible that this is due to the cross-correlation func- 
tion estimator having larger errors than predicted 
by the Poisson estimate (Landy and Szalay 1993). 
Also, at faint magnitudes in the F855 field, the 
early-late cross-correlation function is comparable 
to the expected variations in the galaxy number 
counts introduced by errors in the dust extinction 
estimates. 



To determine the spatial correlation function, 
an estimate of the redshift distribution is required. 
Estimates of the redshift distribution derived from 
luminosity functions, /c— corrections and evolution 
models are subject to uncertainties as a range 
of models can reproduce the observed number 
counts. Models which assume a shape for the red- 
shift distribution (Baugh and Efstathiou 1993) 
are useful for single band imaging of low redshift 
galaxies but are not as effective for multicolor se- 
lected samples at higher redshifts where the selec- 
tion criteria and A;— corrections skew the redshift 
distribution. 

Photometric redshifts contain information on 
the redshift distribution but assume that galax- 
ies with the same colors and magnitudes are at 
the same redshifts. If the errors of the photomet- 
ric redshifts are dominated by the redshift distri- 
bution of galaxies with the same multicolor pho- 
tometry, it should be possible to derive a redshift 
distribution using the measured errors of the pho- 
tometric rcdshifts. The redshift distribution has 
been estimated by smoothing the photometric red- 
shift distribution with a Gaussian where the a of 
the Gaussian is given by the rms of the errors of 
the redshift estimates. At z < 0.05 the redshift 
distribution has been multiplied by 2/0.05 to pre- 
vent an infinite density of galaxies at z = 0. Fig- 
ure 6 shows the observed redshift distribution of 
galaxies from Glazebrook et al. (1995) detected in 
the stacked scans and models derived from their 
photometric redshifts. There is reasonable agree- 
ment between the measured and model distribu- 
tions, though the galaxy number counts are lim- 
ited. While the exact redshift distributions have 
not been determined, it is unlikely that large er- 
rors dominate the model redshift distributions. 

The estimated redshift distributions of the 
early and late subsamples are plotted in Fig- 
ure 7. The photometric redshifts of blue galax- 
ies have large error estimates and this results 
in the late subsample redshift distribution be- 
ing significantly broader than the early subsample 
redshift distribution. The width of the distribu- 
tion is consistent with late-type galaxies having 
smaller fc— corrections and a higher fraction of 
dwarf galaxies than early- type galaxies. It should 
be noted that while the photometric redshifts are 
complete to Bj = 23.5 for the early subsample, 
photometric redshifts for Bj > 22.5 late subsam- 
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pic galaxies arc incomplete and this may slightly 
skew the redshift distribution. 

The amplitudes of the early and late autocor- 
relation functions and early-late cross-correlation 
function have been plotted in Figure 8. To al- 
low the comparison of the amplitude at different 
magnitudes, 7 has been fixed for all magnitudes 
to the average value of 7 for the SGP and F855 
18.0 < Bj < 23.5 correlation functions. Models 
with the clustering fixed in comoving coordinates 
have been fitted to the 18.0 < Bj < 22.5 mag- 
nitude range for the SGP and F855 fields. The 
models have been fitted to this magnitude range as 
it includes a large number of galaxy-galaxy pairs 
while avoiding systematic errors which could be 
present near the magnitude limits of the data. Ad- 
ditional estimates of ro with clustering fixed in co- 
moving and physical coordinates are provided in 
Table 5. 

Figure 8 shows the clustering in the SGP is con- 
sistently stronger than the clustering in F855 for 
18.0 < Bj < 22.5 galaxies. Before concluding 
that large-scale-structure is responsible for the ob- 
served difference between the clustering in the two 
fields, it is useful to exclude the most plausible 
sources of error. Systematic errors in the Bj — R 
colors of ~ 0.2 magnitudes in the two fields only 
slightly alter the observed clustering in each field 
and are inconsistent with galaxy number counts, 
color-color diagrams and published photometry. 
Large uniform errors in the photometry in all 4 
bands can bring the clustering at bright magni- 
tudes into agreement but are also inconsistent with 
galaxy number counts and published photometry. 
The Schlegel et al. (1998) dust maps estimate 
E{B-V) ~ 0.05 in the F855 field and errors in the 
dust extinction estimates could introduce system- 
atic errors. However, errors in the dust extinction 
estimates should produce the largest discrepencies 
at faint magnitudes where the amplitude of the an- 
gular correlation function is small rather than at 
bright magnitudes. 

It is therefore not implausible that the differ- 
ence between the two fields is due to large-scale- 
structure. The distribution of galaxy clusters in 
both fields shows evidence of structures compara- 
ble to the field-of-view (Brown et al. 2000; Brown 
2000) and the SGP may contain several "sheets" of 
galaxies (Broadhurst et al. 1990) though it is un- 
certain if this is particularly unusual (Kaiser and 



Peacock 1991). The effects of large-scale-structure 
on estimates of the correlation function should de- 
crease with increasing survey volume and this is 
consistent with the observed difference between 
the clustering in the two fields decreasing with 
increasing depth. While it is unexpected to see 
weakly clustered late subsample galaxies showing 
the effects of large-scale-structure, the discrepancy 
between the two fields disappears for bluer galax- 
ies selected with U-Bj < 0.4 (Brown et al. 2000). 
As it is uncertain which field is more representa- 
tive of the Universe, the discussion of AGN-galaxy 
clustering in the SGP assumes that estimates of 
ro could have systematic errors comparable to the 
difference between the estimates of ro in the SGP 
and F855 (- Ihr^M^c). 

Models for the clustering of early subsample 
galaxies in each field are a reasonable approxima- 
tion to the observed clustering and the estimates 
of the spatial correlation function are compara- 
ble to the clustering of early-type galaxies in the 
local Universe (Loveday et al. 1995; Guzzo et al. 
1997). In contrast, the observed amplitude of the 
late subsample autocorrelation function is an or- 
der of magnitude weaker and the models arc a 
poor fit to the data. At Bj > 22 the ampli- 
tude of the correlation function rapidly decreases 
due to the increasing number of weakly clustered 
faint blue galaxies (Efstathiou et al. 1991; Brown 
et al. 2000). The amplitude of the early-late 
cross-correlation function also decreases at faint 
magnitudes though models of the spatial corre- 
lation function are better approximations to the 
observed clustering. While the values of ro for the 
early-late cross-correlation function are similar to 
those for the late subsample autocorrelation func- 
tion, the larger value of 7 results in significantly 
stronger clustering on scales of < l/i~^Mpc. 

5. AGN-Galaxy Clustering 

Estimates of the cross-correlation of AGN with 
the early and late subsamplcs are summarized in 
Table 6. The AGN-early cross-correlation is sig- 
nificantly stronger on scales > 1' while on smaller 
scales the signal-to-noise is poor. Despite the 
strength of the AGN-early cross-correlation, in- 
spection of Table 6 shows most AGN-galaxy pairs 
are AGN-late pairs and the cross-correlation of 
AGN with all galaxies will be similar to the AGN- 
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late cross-correlation function. 

The AGN-early and AGN-late angular cross- 
correlation functions for Bj < 22.5 galaxies are 
shown in Figure 9. Power-laws with 7 fixed to 
1.90 and 1.65 have been fitted to the AGN-early 
and AGN-late data on angular scales < 0.3°. The 
AGN-early cross-correlation function is strong and 
the power-law fit is a good approximation to the 
observed clustering though most of the data points 
are > la from the fit. In contrast, the AGN-late 
cross-correlation is weak and a power-law is a poor 
fit to the data. From the available data, the am- 
plitude of the AGN-late cross-correlation function 
can not be reliably measured. 

The amplitude of the AGN-early cross-correlation 
function is plotted in Figure 10 along with mod- 
els of the cross-correlation with evolution fixed in 
comoving and physical coordinates. While the am- 
plitude of the angular cross-correlation function is 
typically less than the early subsample autocorre- 
lation function, the difference between the early 
subsample and AGN redshift distributions results 
in the estimates of tq being significantly higher 
than that of the early subsample autocorrelation 
function. 

The errors of the lli{0) estimates are dominated 
by galaxies along the line-of-sight which are not 
associated with the AGN. The signal and sig- 
nificance of the correlation function can be im- 
proved if the number of unassociated galaxies in- 
cluded in the estimate of u){9) are reduced. Early- 
type galaxies have comparatively small photomet- 
ric redshift errors so it is possible to define a nar- 
row range of photometric redshifts which could 
be associated with the AGN. Figure 11 shows the 
AGN-early cross-correlation function for galaxies 
with photometric redshifts within 2a of the AGN 
redshift. The inclusion of photometric redshift 
constraints has significantly improved the signal- 
to-noise and the observed correlation function is 
more consistent with a power-law than Figure 9. 
While deviations from a power-law have been seen 
in some estimates of AGN-galaxy correlation func- 
tions (Hall and Green 1998: Groom and Shanks 
1999), there are no statistically significant breaks 
from a power-law on scales < 0.3° in Figure 11. 

The spatial AGN-early cross-correlation can bo 
determined by deriving the redshift distribution of 
galaxies with photometric redshifts within 2a of 
each AGN when applying Limber's equation. The 



amplitude of the AGN-early cross-correlation and 
models of the spatial cross-correlation are shown 
in Figure 12 and the fit to the Bj = 22.5 data 
point is weaker than in Figure 10. However, the 
estimate of ro and the power-law form of the angu- 
lar correlation function in Figure 1 1 are consistent 
with AGN environments being similar to the envi- 
ronments of strongly clustered galaxies such as the 
early subsample and low redshift ellipticals (Guzzo 
et al. 1997). 

The estimates of ro are dependent on the as- 
sumed value of 7 and it is possible that 7 ^ 
1.90. If radio-quiet AGN were in spiral galax- 
ies (Hutchings et al. 1984; Malkan et al. 1984) 
and had similar clustering properties to late sub- 
sample galaxies, the value of 7 for the AGN-early 
cross-correlation function would be ~ 1.65. How- 
ever, when the AGN-early cross-correlation func- 
tion was fitted with power-laws with 7 < 1.65, the 
estimates of ro increased to > ll/i~^Mpc which 
is approximately double the estimate of ro for the 
early-late spatial cross-correlation function. 

There may be a correlation between AGN lu- 
minosity and host galaxy if luminous AGN oc- 
cur more frequently in early-type galaxies than 
Seyferts (McLure et al. 1999). If this is the case, 
the correlation between galaxy morphology and 
environment may result in a correlation between 
AGN luminosity and environment. Figure 13 com- 
pares the amplitudes of the Seyfert 1 {Mb, < 
-21.5) and QSO {Mb,, > -21.5) cross-correlation 
functions with the early subsample. Fits to the 
Bj = 22.5 data point show Seyfert Is are corre- 
lated with slightly richer environments but differ- 
ence between the ro estimates has < la signifi- 
cance. 

6. Discussion 

While the detection of radio-quiet AGN in en- 
vironments similar to early-type galaxies may ap- 
pear to contradict previous studies at similar red- 
shifts (Yec and Green 1987; Boyle and Couch 
1993; Groom and Shanks 1999; Teplitz et al. 1999; 
Smith et al. 2000), this assumes all studies of 
AGN-galaxy clustering have used the same galaxy 
population to measure AGN environments. How- 
ever, as color selection has been used to select 
galaxies for this study, the properties of the early 
and late subsamples will differ significantly from 
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catalogues of galaxies which are selected by ob- 
served broadband flux. 

Early subsample galaxies contribute ~ 30% of 
the total galaxy counts at Bj < 21 and only 
~ 15% of the counts at Bj > 22. If the AGN-late 
cross-correlation function is assumed to be neg- 
ligible and the AGN-early cross-correlation func- 
tion is a power law with ro ^ 9/i~^Mpc, then the 
AGN-galaxy cross-correlation function should be 
a power-law with rg ^ 5ft-"^Mpc at Bj < 21 and 
ro ~ 3/i~^Mpc at Bj > 22. In redder bands where 
the early-type galaxies comprise a larger fraction 
of the total galaxy number counts, stronger clus- 
tering will be measured. Early subsample galax- 
ies comprise ~ 25% of the total R < 21.5 galaxy 
counts and the AGN-galaxy correlation function 
would be expected to be ro ~ 4/i~^Mpc even 
though the redshift range sampled is compara- 
ble to that of Bj ~ 23 galaxies. These esti- 
mates of ro for Bj and R limited samples are 
consistent with measurements of the AGN-galaxy 
cross-correlation function derived from single band 
imaging by Yee and Green (1987), EUingson et al. 
(1991) and Smith et al. (1995, 2000). 

While radio-quiet AGN occur in environments 
comparable to early- type galaxies, this should not 
be interpreted as a direct correlation between rich- 
ness and AGN activity. As shown in Figure 13, 
there is no evidence for a strong correlation be- 
tween AGN luminosity and environment. HST 
imaging of radio-quiet QSOs and X-ray selected 
AGN indicates ^ 75% of host galaxies have mor- 
phologies earlier than Sbc galaxies (McLure et al. 
1999; Schade et al. 2000). It is therefore plausi- 
ble that AGN activity is not significantly effected 
by the host environment and that the observed 
strength of the AGN-early cross-correlation func- 
tion is due to the correlation between galaxy mor- 
phology and clustering properties. 

The strength AGN-early cross-correlation func- 
tion implies that it should be detectable at z ^ 
0.5. Using samples of early-type galaxies obtained 
with CCD mosaics on 4-m class telescopes, it will 
be possible to measure the evolution of the AGN- 
early cross-correlation function at z > 1. This will 
provide an estimate of the bias of AGN with re- 
spect to galaxies which are used to measure large- 
scale-structure at lower redshifts. If the cluster- 
ing properties of AGN are due to the distribu- 
tion of host galaxy morphologies, the evolution of 



the AGN-early cross-correlation function will also 
place strong constraints on the properties of AGN 
host galaxies as a function of redshift. 

7. Summary 

The environments of galaxies and 69 0.2 < z < 
0.7 AGN have been measured using photometric 
redshifts and color criteria to select galaxy types 
to Bj ~ 23.5. The key conclusions are 

(i) The clustering of early subsample galaxies is 
strong across the observed magnitude range with 
ro ~ 7h-^Mpc and 7 ~ 1.90. 

(ii) The autocorrelation function of late subsample 
galaxies is weak and decreases to ro < 4/i~-^Mpc 
at i? J ~ 23.5. This is probably due to the increas- 
ing fraction of weakly clustered blue galaxies at 
faint apparent magnitudes. 

(iii) The cross-correlation function of radio-quiet 
-24 < Mb < -19 AGN with early subsample 
galaxies has been detected with high significance 
on scales < 1°. The AGN-early spatial cross- 
correlation function is stronger than the early 
subsample autocorrelation function and is com- 
parable to the clustering of elliptical galaxies at 
2; ~ 0. 

(iv) The AGN-late cross-correlation function is 
very weak and has been detected with low sig- 
nificance. As the fraction of late-type galaxies in 
magnitude limited samples increases with survey 
depth, the cross-correlation function of AGN with 
all galaxies decreases with increasing magnitude. 

(v) The signal-to-noise of AGN-galaxy cross- 
correlation functions is significantly improved 
by using photometric redshifts to reject galaxies 
which can not be associated with the AGN from 
the correlation function estimate. 

(vi) The correlation between AGN optical lumi- 
nosity and host environment is weak and has not 
been detected at a significant level in this work. 
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Fig. 1. — The absolute magnitudes and redshifts 
of the 0.2 < z < 0.7 AGN selected from the 
La Franca ct al. (1999) catalogue. The absolute 
magnitudes have been determined with fl = 1, 
Hq = 100kms~^Mpc~^ and fc— corrections approx- 
imated by k{z) = —0.42;. 




J I I 1 I I I 1 L_ 

0.2 0.4 0.6 0.8 



Z — spectroscopic 



Fig. 2. — Comparison of photometric and spectro- 
scopic redshifts of galaxies. The accuracy of the 
redshift estimates derived with 3 or more colors 
is significantly better than redshifts derived with 



This 2-column preprint was prepared with the AAS lAT^X BjR photometry, 

macros v5.0. 



9 




Spectroscopic Redshift 



Fig. 3. — Comparison of photometric and spec- 
troscopic rcdshifts for galaxies with colors redder 
than non-evolving Sbc galaxies. The dashed lines 
show the ±1<T error estimates of the photometric 
redshifts. While some of the galaxies arc only de- 
tected in Bj and R, the correlation between pho- 
tometric and spectroscopic redshift is still good. 



Fig. 4. — The colors and photometric redshifts 
of galaxies with morphological classifications from 
Abraham et al. (1996) and Small ct al. (1997). 
The curved line is an estimate of the color of a 
non-evolving Sbc determined with a polynomial fit 
to the fc— corrections from Coleman ct al. (1980). 
Most early-type galaxies are redder than the non- 
evolving Sbc. 
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F855 Early Power Law Fit •y=1.96 




J SGP Late Data with 1(7 Error Bars 

f FB55 Late Data with Iff Error Bars 

— SGP Late Power Law Fit 7=1.32 

— F855 Late Power Law Fit 7=1.38 
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Fig. 6. — A comparison of the observed and model 
redshift distributions of galaxies from Glazebrook 
et al. (1995) which were detected in the stacked 
scans. The left panel shows all galaxy types while 
the right panel shows galaxies with colors redder 
than Sbc galaxies. While limited by the small 
number of galaxies available, the models are simi- 
lar to the observed redshift distributions with both 
peaking at similar redshifts. 



0.1 

6 (Degrees) 



Fig. 5.— The 18.0 < Bj < 22.5 galaxy angu- 
lar correlation functions for the SGP and F855 
fields. The autocorrelation function of early sub- 
sample galaxies is significantly stronger than the 
autocorrelation function of late subsample galax- 
ies. Power laws fitted to the < 0.3° data are good 
approximations to the observed clustering on most 
angular scales. 
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Fig. 7. The estimated redshift distribution of 
the early and late subsamples in the SGP. The 
curves range from 18.0 < Bj < 20.0 to 18.0 < 
B,j < 23.5 in half magnitude steps. The large 
redshift errors for blue galaxies results in the late 
subsample redshift distribution being significantly 
broader than the early subsample redshift distri- 
bution. At magnitudes fainter than Bj = 22.5 
photometric redshifts are not available for all late 
subsample galaxies and this may produce errors in 
the model redshift distribution. 




E 



n-0.3. A-0.7 
] SGP Early Data with Ict Error Bars 
^ F855 Early Data with 1(j Error Bars 

SGP Early Model r„=6.6 h~' Mpc, 7'=1.90, £=-1.10 

- F855 Early Model r^=5.3 h"' Mpc, y=T.90, e=-l.10 



n-0.3, A-0.7 • , 

f SGP Late Data with Icr Error Bars 
J F855 Late Data with 1(j Error Bore 

SGP Late Model r„=4.5 h"' Mpc, 7=1.35, e=-1.65 

— FB55 Late Model r^=3.5 h~' Mpc. 7=1.35, e=-1.65 



E 

o d 




n-0.3, A-0.7 
I SGP Early-Late Data with la Error Bars 
f F855 Early-Late Data with la Error Bars 

SGP Early-Late Model r„=5.0 h~' Mpc, 7=1.65, e=-1.35 : 

— F855 Early-Late Model r|,=3.8 h~' Mpc, 7=1.65, e=-1.35' 



Fig. 8. — The amplitudes of the early and late 
autocorrelation functions and the early-late cross- 
correlation function. Models fitted to the Bj = 
22.5 data points with clustering fixed in comoving 
coordinates are shown. The median redshift of the 
sample as a function of depth is listed at the top 
of the plot. The difference between the observed 
clustering in the two fields generally decreases as 
a function of magnitude. 
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Fig. 9 — The AGN-early and AGN-late angular 
cross-correlation functions. The AGN-early cor- 
relation function is significantly strongc^r than the; 
AGN-late correlation function. While a power-law 
is a good approximation to the observed AGN- 
early cross-correlation function, a power-law fit to 
the AGN-late cross-correlation function is a poor 
fit to the data. 
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Fig. 10. — The amplitude of the AGN-early angu- 
lar cross-correlation function. Models with clus- 
tering fixed in comoving and physical coordinates 
have been fitted to the Bj = 22.5 data point. The 
clustering is significantly stronger than the auto- 
correlation function of early subsample galaxies. 
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Fig. 11. — The angular cross-correlation of AGN 
and Bj < 22.5 early subsample galaxies with pho- 
tometric redshift constraints applied to the pair 
counts. Poisson statistics have been used to de- 
termine the la errors shown with the data points. 
A power law with 7 = 1.90 is a good fit to the 
data. 
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Fig. 12.— The amplitude of the AGN-early 
cross-correlation function with photometric red- 
shift constraints applied. The errors arc signifi- 
cantly reduced and the best fits to the data are 
slightly decreased. However, the amplitude of the 
clustering is consistent with AGN being in envi- 
ronments comparable to elliptical galaxies. The 
estimated median redshifts are for early subsam- 
plc galaxies that satisfy the photometric rcdshift 
constraints rather than the entire early subsample. 
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Fig. 13. — The amplitude of the Seyfert 1 and 
QSO cross-correlation functions with photometric 
redshift constraints. There is no evidence for a 
strong correlation between AGN luminosity and 
host environment. 
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Table 1 

a sample of studies of agn environments 



Reference 


AGN 


Number 


AGN 


AGN 




Galaxy 


Correlation 




Type 


of 


Redshift 


Magnitudes 


Magnitude 


Function 






AGN 


Range 






Limit 


Estimate 


De Robertis et al. (1998) 


Seyferts 


33 


z < 0.04 


-23.0 < Mr < 


-18.5 


Mr < -19.5 


ro ~ 5.4?i~'^Mpc 


Laurikaninen and Salo (1995) Seyfert 1 


55 


z < 0.05 


-22.0 < Mv < 


-17.0 


Bj < 21 


ro ~ 5/i~^Mpc 


Laurikaninen and Salo (1995) Seyfert 2 


49 


z < 0.05 


-21.5 < Mv < 


—21.5 


Bj < 21 


ro > 5/i~^Mpc 


Yee and Green (1987) 


RQ QSOs 


33 


z < 0.3 


-26.0 <Mr< 


-23.5 


r < 21.5 


ro ~ 8/i~^Mpc 


Smith et al. (1995) 


X-ray QSOs 


169 


z < 0.3 


-23.5 < Mv < 


-17.0 


Bj < 20.5 


7*0 ^ 5/i~^Mpc 


Yee and Green (1987) 


RQ QSOs 


7 


0.3 < z < 0.6 


-26.5 < Mr < 


-23.5 


r < 21.5 


ro ~ 5/i~^Mpc 


Ellingson et al. (1991) 


RQ QSOs 


46 


0.3 < z < 0.6 


-24.0 < Mr < 


-18.5 


r 23.5 


ro ~ 6/(~^Mpc 


This work (AGN-Early) 


UBR AGN 


69 


0.2 < z < 0.7 


-24.0 < Mb < 


-19.0 


Bj < 23.5 


ro ~ 8/i~^Mpc 


Smith et al. (2000) 


X-ray QSOs 


83 


0.3 < z < 0.7 


-26.0 < Mv < 


-17.5 


V = 23 


ro ~ 3/(~^Mpc 


Boyle and Couch (1993) 


QSOs 


27 


0.9 < z < 1.7 


-24.0 < Mb < 


-22.0 


R <23 


No correlation 


Tcplitz ct al. (1999) 


RQ QSOs 


30 


0.9 < z < 2.1 


-26.5 < Mv < 


-20.6 


I,J,Hr^ 21, K 


2(T correlation 


Groom and Shanks (1999) 


RQ QSOs 


150 


0.0 < z < 3.2 


B <22 




Bj < 23 


No correlation 


Yee and Green (1987) 


RL QSOs 


10 


0.3 < z < 0.5 


-25.0 < Mr < 


-23.5 


r < 21.5 


ro ~ 9/i~^Mpc 


Yee and Green (1987) 


RL QSOs 


9 


0.55 < z < 0.65 


-25.5 < Mr < 


-23.5 


r < 21.5 


ro ~ 17?i~^Mpc 


Wold cl, al (2000) 


RL QSOs 


21 


n.5n < z < 0,82 


-24, .'i < Mb < 


-19.!) 


V. n ~ 2.'!. 5. / 


r„ ~ ll,7/)"iMpr 


Hall and C4r<-<-u [I'.m^) 


liL (JSOs 


:il 


1,(1 < < 2,(1 


-27, -i . . ..\U < 


-21.(1 


A > 1') 


liicli, (1 cliL.Mcr,^ 



Table 2 

The Early Subsample Angular Correlation Function 



Field 


SGP 






F855 




Magnitude Range 


Ngal 7 


LJ{1') X 10^ 


Ngal 


7 


X 10^ 




2206 1.70 ±0.11 


1352 ± 233 


2599 


1.79 ±0.31 


782 ± 196 


18.0 <Bj< 21.0 


5187 1.79 ±0.08 


830 ± 83 


5553 


2.02 ±0.19 


525 ± 78 


18.0 <Bj < 22.0 


11056 1.74 ±0.04 


510 ±51 


11303 


1.94 ±0.09 


326 ± 34 


18.0 <Bj < 23.0 


21833 1.81 ±0.05 


283 ± 28 


21652 


1.87 ±0.07 


194 ± 19 


18.0 < Bj < 23.5 


29290 1.84 ±0.06 


221 ± 22 


27915 


1.89 ±0.05 


169 ± 17 






Table 3 








The Late Subsample Angular Correlation Function 


Field 


SGP 






F855 




Magnitude Range 


Ngal 7 


w(l') X 10^ 


Ngal 


7 


W(l') X 103 


i8n< R7<2nn 


3854 1.35 ±0.27 


443 ± 257 


4648 


1.57 ±0.37 


289 ± 87 


18.0 < Bj < 21.0 


11952 1.39 ±0.13 


269 ± 25 


13866 


1.62 ±0.16 


157 ±21 


18.0 <Bj< 22.0 


37785 1.33 ±0.10 


120 ±9 


43047 


1.42 ±0.15 


76 ±7 


18.0 <Bj< 23.0 


128278 1.29 ±0.09 


44 ±8 


134331 


1.31 ±0.18 


37 ±16 


18.0 < Bj < 23.5 


207390 1.34 ±0.12 


28 ±4 


183347 


1.33 ±0.12 


27 ±5 






Table 4 








The Early-Late Angular Cross-Correlation Function 


Field 


SGP 




F855 






Magnitude Range 


7 toil') X 10^ 7 


L0{1') X 103 




18.0 <Bj < 20.0 


1.43 ±0.25 430 ± 105 1.33 ±0.35 ; 


202 ± 70 




18.0 <Bj < 21.0 


1.67 ±0.08 322 ±30 1.48 ±0.14 ; 


204 ± 22 




18.0 < Bj < 22.0 


1.73 ±0.06 183 ±10 1.58 ±0.08 


102 ±8 




18.0 <Bj< 23.0 


1.76 ±0.04 85 ±4 1.59 ±0.08 


52 ±3 




18.0 <Bj < 23.5 


1.65 ±0.04 56 ±3 1.46 ±0.09 


40 ±2 
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Table 5 

The spatial autocorrelation and cross-correlation functions of galaxies. 



Correlation 


Galaxy Magnitude 


Field 


7 


ro (/i~^Mpc) 


ro (/i~^Mpc) 


Function Samples 


Range 






e = 1 — 7 


e = 


Early-Early 


18.0 <Bj < 22.5 


SGP 


1.90 


6.6 ±0.3 


7.7 ±0.3 


Early-Early 


18.0 < Bj < 22.5 


F855 


1.90 


5.3 ±0.3 


6.1 ±0.3 


Late-Late 


18.0 <Bj < 22.5 


SGP 


1.35 


4.5 ±0.2 


5.9 ±0.3 


Late-Late 


18.0 <Bj < 22.5 


F855 


1.35 


3.5 ±0.2 


4.7 ±0.3 


Early-Late 


18.0 < Bj < 22.5 


SGP 


1.65 


5.0 ±0.2 


6.1 ±0.2 


Early-Late 


18.0 <Bj < 22.5 


F855 


1.65 


3.8 ±0.2 


4.7 ±0.2 


Early-Early 


18.0 < Bj < 23.5 


SGP 


1.90 


5.8 ±0.2 


7.0 ±0.2 


Early-Early 


18.0 <Bj < 23.5 


F855 


1.90 


5.2 ±0.3 


6.2 ±0.3 


Late-Late 


18.0 <Bj < 23.5 


SGP 


1.35 


2.5 ±0.2 


3.6 ±0.3 


Late- Late 


18.0 <Bj < 23.5 


F855 


1.35 


2.5 ±0.2 


3.5 ±0.3 


Early-Late 


18.0 <Bj < 23.5 


SGP 


1.65 


4.0 ±0.1 


5.1 ±0.2 


Early-Late 


18.0 <Bj < 23.5 


F855 


1.65 


3.0 ±0.2 


4.0 ±0.2 



Table 6 

The cross-correlation function of 0.2 < z < 0.7 AGN and early subsample galaxies. The 
values of u{9) have been corrected for the integral constraint which is ~ 10% OF the 

ESTIMATE OF THE CORRELATION FUNCTION AT 1'. 



AGN-Early 



Magnitude Range 


Anj 


gle 


AG 


AR X 20 




18.0 < Bj 


< 22.5 


10" <e < 20" 


6 


81 


0.48 ± 0.84 


18.0 < Bj 


< 22.5 


20" < 6 


' < 40" 


15 


281 


0.07 ±0.39 


18.0 < Bj 


< 22.5 


40" < 6 


' < 80" 


59 


1026 


0.15 ±0.15 


18.0 < Bj 


< 22.5 


80" < 6 


< 160" 


272 


4299 


0.27 ±0.07 


18.0 < Bj 


< 22.5 


160" < 6 


' < 320" 


882 


17203 


0.03 ±0.03 


18.0 < Bj 


< 22.5 


320" < 6 


' < 640" 


3704 


68428 


0.08 ±0.02 


18.0 < Bj 


< 23.5 


10" < 6 


' < 20" 


10 


132 


0.52 ±0.61 


18.0 < Bj 


< 23.5 


20" < 6 


' < 40" 


28 


520 


0.08 ±0.21 


18.0 < Bj 


< 23.5 


40" <e < 80" 


114 


1956 


0.17±0.11 


18.0 < Bj 


< 23.5 


80" <e < 160" 


497 


8128 


0.04 ±0.06 


18.0 < Bj 


< 23.5 


160" < 6 


' < 320" 


1676 


32371 


0.07 ±0.03 


18.0 < Bj 


< 23.5 


320" <e < 640" 


6881 


129220 


0.02 ±0.01 
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Table 7 

The cross-correlation function of 0.2 < ^; < 0.7 AGN and late subsample galaxies. 



AGN-Late 

Magnitude Range Angle AG AR x 20 Ld{6) 



18.0 < 


Bj 


< 


22.5 


10" < 1 


9 < 20" 


15 


320 


-0.06 ±0.35 


18.0 < 


Bj 


< 


22.5 


20" < 1 


9 < 40" 


53 


1183 


-0.10 ±0.12 


18.0 < 


Bj 


< 


22.5 


40" < 1 


9 < 80" 


263 


4680 


0.13 ±0.07 


18.0 < 


Bj 


< 


22.5 


80" < 6 


1 < 160" 


949 


19111 


0.00 ±0.03 


18.0 < 


Bj 


< 


22.5 


160" < 1 


9 < 320" 


3812 


76187 


0.00 ±0.02 


18.0 < 


Bj 


< 


22.5 


320" < 1 


9 < 640" 


15602 


303404 


0.03 ±0.01 


18.0 < 


Bj 


< 


23.5 


10" < ( 


9 < 20" 


42 


905 


-0.07 ±0.14 


18.0 < 


Bj 


< 


23.5 


20" < ( 


9 < 40" 


175 


3333 


0.05 ± 0.08 


18.0 < 


Bj 


< 


23.5 


40" < 1 


9 < 80" 


751 


13747 


0.09 ±0.04 


18.0 < 


Bj 


< 


23.5 


80" < 6 


• < 160" 


2805 


55104 


0.02 ± 0.02 


18.0 < 


Bj 


< 


23.5 


160" <e< 320" 


11166 


221132 


0.01 ±0.01 


18.0 < 


Bj 


< 


23.5 


320" <e < 640" 


45032 


880067 


0.03 ±0.01 



Table 8 

The cross-correlation of AGN with the early subsample with photometric redshift 

rONSTRAIXTR APPLIED TO THE DATA. 



Magnitude Range 


Angle 


AG 


AR X 20 


^{0) 


18.0 < Bj 


< 22.5 


10" <e < 20" 


1 


20 


0.07 ±2.20 


18.0 < Bj 


< 22.5 


20" <e < 40" 


10 


133 


0.51 ±0.61 


18.0 < Bj 


< 22.5 


40" <e < 80" 


38 


518 


0.48 ±0.23 


18.0 < Bj 


< 22.5 


80" <9 < 160" 


126 


2104 


0.21 ±0.10 


18.0 < Bj 


< 22.5 


160" <e < 320" 


438 


7947 


0.11 ±0.05 


18.0 < Bj 


< 22.5 


320" <d < 640" 


1704 


31708 


0.08 ±0.03 


18.0 < Bj 


< 23.5 


10" <e < 20" 


5 


71 


0.41 ±0.91 


18.0 < Bj 


< 23.5 


20" <0 < 40" 


23 


373 


0.24 ±0.26 


18.0 < Bj 


< 23.5 


40" <e < 80" 


89 


1388 


0.28 ±0.13 


18.0 < Bj 


< 23.5 


80" <9 < 160" 


340 


5637 


0.21 ±0.07 


18.0 < Bj 


< 23.5 


160" < 6* < 320" 


1172 


21818 


0.08 ±0.03 


18.0 < Bj 


< 23.5 


320" <e< 640" 


4661 


87823 


0.08 ±0.02 
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